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1. Introduction 

The Casimir effect [lj can be seen as an effect of the zero point energy of vacuum 
which emerges due to the non-commutativity of quantum operators upon quantisation 
of the electromagnetic (EM) field. Although formally infinite in magnitude, the EM 
field density in bulk undergoes finite alterations when dielectric or metal boundaries 
are introduced in the system, giving rise to finite and measurable forces. As is well 
known, at nanometre to micrometre separations the Casimir attraction between bodies 
becomes significant, and the effect has attracted much attention during the last decade 
in the wake of the rapid advances in nanotechnology. The existence of the Casimir force 
was shown experimentally as early as 1958 by Spaarnay [2], yet only recently new and 
much more precise measurements of Lamoreaux and others (see the review [3]) have 
boosted the interest in the effect from a much broader audience. Experiments like that 
of Mohideen and Roy [I], and the very recent one of Harber et. al. [5], making use 
of the oscillations of a magnetically trapped Bose Einstein condensate, were subject 
to widespread regard. The same was true for the nonlinear micromechanical Casimir 
oscillator experiment of Chan et al. [6J Ej . 

Recent reviews on the Casimir effect are given in Refs. [3], El El OH EI]- Much 
information about recent developments can also be found in the special issues of J. 
Phys. A: Math. Gen. (May 2006) p2], and of New J. Phys. (October 2006) p3]. 

Actual calculations of Casimir forces are usually performed via two different routes 
[8]; either by summation of the energy of discrete quantum modes of the EM field (cf., 
for instance, [H]), or via a Green's function method first developed by Lifshitz [15]. 
Mode summation, despite its advantage of a simpler and more transparent formalism, 
is usually far inferior. In practice it is only in systems where quantum energy states are 
known that energy summation can be carried out explicitly. This requires the system to 
be highly symmetric, and favour assumptions such as perfectly conducting walls like in 
the original Casimir problem. Geometries in which quantum states are known exactly, 
unfortunately, are few. 

The method of calculating the force through Green's functions avoids some but 
not all of these problems; exact solutions are still only known in highly symmetrical 
systems such as infinitely large parallel plates or concentric spheres. Via the fluctuation- 
dissipation theorem, the EM field energy density is linked directly to the photonic 
Green's function, and the force surface density acting on boundaries can be calculated, 
at least in principle. The theory of Green's functions and the application of them will 
be central in the present paper. 

The purpose of the present work is twofold. First, we intend to explore some of the 
delicate issues that occur in the Green's function formalism in typical settings involving 
dielectric boundaries. Upon relating the two-point functions to the Green's function one 
may choose to calculate the Green function in full [SI [16] . The method is complete but 
may appear cumbersome, at least so in the presence of several dielectric surfaces. It is 
possible to reduce the calculational burden somewhat by simplifying the Green function 
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expressions, by omitting those parts that do not contribute to the Casimir force. This 
means that one works with "effective" Green functions. This method is employed and 
briefly discussed by Lifshitz and co-workers, cf. e.g. [IT]. The connections between the 
different kinds of Green's functions are in our opinion far from trivial, and we therefore 
believe it of interest to present some of the formulas that we have compiled and which 
have turned out to be useful in practice. 

As for the calculational technique for the Casimir force in a multilayer system, 
there exists a powerful formalism worked out, in particular, by Tomas [18]. In turn, 
this formalism was based on work by Mills and Maradudin two decades earlier [19J. 
One of us recently made a review of this technique, with various applications . We 
shall make use of this technique in the following. In company with the by now classic 
theory of Lifshitz and co-workers [T5l [21] and the standard Fresnel theory in optics, the 
necessary set of tools is provided. 

Our second purpose is to apply the formalism to concrete calculations of the Casimir 
pressure on a dielectric plate in a multilayer setting. Especially, we will consider the 
pressure on a plate situated in a cavity (5-zone system). We work out force expressions 
and eigenfrequency changes when the plate is acted upon by a harmonic-oscillator 
mechanical force (spring constant k) in addition to the Casimir force, and is brought to 
oscillate horizontally. To our knowledge, explicit calculations of this sort have not been 
made before. A chief motivation for this kind of calculation is that we wish to evaluate 
the magnitudes of thermal corrections to the Casimir pressure. In recent years there 
have been lively discussions in the literature about the thermal corrections; for some 
statements of both sides of the controversy, see Refs. [22l [23l [2^1 125| f26| f27j f28l [29] . 
We hope that the consideration of planar multilayer systems may provide additional 
insight into the temperature problem. 

We will be considering uniformly heated systems only. The recent experiment of 
Harber et al. [5] investigated the surface-atom force at thermal equilibrium at room 
temperature, the goal being to measure the surface-atom force at very large distances, 
taking into account the peculiar properties of a Bose-Einstein condensate gas. Later, 
the same group investigated the non-equilibrium effect [30] . This paper seems to have 
reported the first accurate measurement of the thermal effect (of any kind) of the Casimir 
force, in good agreement with earlier theoretical predictions [31] (cf. also the prior theory 
of Pitaevskii on the non-equilibrium dynamics of EM fluctuations [32]). Consideration 
of such systems lies, however, outside the scope of the present paper. 

The following point ought also to be commented upon, although it is not a chief 
ingredient of the present paper: Our problem bears a relationship to the famous 
Abraham-Minkowski controversy, or more generally the question of how one should 
construct the correct form of the EM energy-momentum tensor in a medium. This 
problem has been discussed more and less intensely ever since Abraham and Minkowski 
proposed their energy-momentum expressions around 1910. The advent of accurate 
experiments, in particular, has aided a better insight into this complicated aspect of 
field-matter interacting systems. Some years ago, one of the present authors wrote a 
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review of the experimental status in the field [33J (cf. also |34J). There is by now a rather 
extensive literature in this field; some papers are listed in [35], [361 EH EH EHl HOj I4TI 142] . 
In the present case, where the EM surface force on a dielectric boundary results from 
integration of the volume force density across the boundary region, the Abraham and 
Minkowski predictions actually become equal. Recently, in a series of papers Raabe and 
Welsch have expressed the opinion that the Abraham- Minkowski theory is inadequate 
and that a different form of the EM energy-momentum tensor has to be employed 
[4"3~1 14*41 1451 146] . We cannot agree with this statement, however. All the experiments 
in optics that we are aware of can be explained in terms of the Abraham-Minkowski 
theory in a straightforward way. One typical example is provided, for instance, by the 
oscillations of a water droplet illuminated by a laser pulse. Some years ago, Zhang and 
Chang made an experiment in which the oscillations of the droplet surface were clearly 
detectable [47] . It was later shown theoretically how the use of the Abraham-Minkowski 
theory could reproduce the observed results to a reasonable accuracy [481 HH] • In our 
theory below, we will use the Abraham-Minkowski theory throughout. 

SI units are used throughout the calculations, and permittivity e and permeability 
// are defined as relative (nondimensional) quantities. We thus write D = eoeE, 
B = /i /iH. 

The outline of the paper is as follows. In the next section we analyse the 5-layered 
magnetodielectric system (figure 1), presenting the full Green's function as well as its 
effective (or reduced) counterpart. We here aim at elucidating some points in the 
formalism that in our opinion are rather delicate. Section 3 is devoted to a study 
of an oscillating slab in a Casimir cavity, permitting, in principle at least, how the 
change in the eigenfrequency of the slab with respect to the temperature can give us 
information about the temperature dependence of the Casimir force. Section 4 discusses 
more extensively the relationships between the Green's two-point functions as introduced 
by Lifshitz et al., and by Schwinger et al. In section 5 we present results of numerical 
Casimir force calculations for selected substances, taking Al as example of a metal, and 
teflon FEP as example of a dielectricl§|. In section 6 we consider the effect of finite slab 
thickness, i.e. the "leakage" of vacuum radiation from one gap to the other. We find the 
striking result that for dielectrics the relative finite thickness correction is much larger 
than for metals. For teflon FEP versus Al the relative correction is almost two orders 
in magnitude higher. 

A word is called for, as regards the permeability fi. As anticipated above, we allow 
fi to be different from 1. This is motivated chiefly by completeness, and is physically 

§ As a word of caution, we mention here that our permittivity data for metals are intended to hold 
in the bulk, whereas in practice the real and imaginary parts of the index of refraction for metals 
come from ellipsometry measurements, and are thus really surface measurements. There is an inherent 
uncertainty in the calculated results coming from this circumstance, of unknown magnitude, although 
in our opinion the corrections will hardly exceed the 1% level due to the general robustness of the force 
expression against permittivity variations. Ideally, information about the permittivity versus imaginary 
frequency would be desirable, for a metallic film. We thank Steve Lamoreaux for comments on this 
point. 
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an idealization. It is known that the permeability for most materials is lossy at high 
frequencies, corresponding to imaginary values for fi. That phenomenon is limited to 
a restricted frequency interval, however, (10 - 100 GHz), and loses effect at the higher 
frequencies. 



2. Casimir force on a slab in a cavity 



We shall consider a 5-layered magnetodielectric system such as depicted in figure UJ The 
analytical calculation of the Casimir force density acting on the slab in such a geometry 
is well known; it may be calculated, quite simply, by a straightforward generalisation 
of the famous calculation by Lifshitz and co-workers used for the simpler, three-layered 
system of two half-spaces separated by a gap [T71 [21] . 

Rather than starting from the photonic Green's function as a propagator as 
known from quantum electrodynamics, we introduce classical and macroscopic two- 
point (Green's) function according to the convention of Schwinger et.al. [52] as 

E(x) = - f dV F (x, x') ■ P(x'), (1) 
eo J 

where x = (r, t). Due to causality, t' is only integrated over the region t' < t. It follows 
from Maxwell's equations that T obeys the relation 

VxVxf (,/;„) - '^p^ ? (r,r» = ^S(r - r>) T, (2) 
where we have performed a Fourier transformation according to 

T(x,x')= ^e— T(r,r';u), (3) 

with r = t — t'. 

A comparison of §2§ with the corresponding equation in [TTI [2Tj shows formally that 
T is essentially equivalent with the retarded photonic Green's function in a mediunjjj]. 
The physical connection is not entirely trivial, however. As motivation we notice that §1§ 
expresses the linear relation between the dipole density at x' and the resulting electric 
field at x, in essence the extent to which an EM field is able to propagate from x' to 
x. This is exactly the classical analogy of the quantum definition of a Green's function 
propagator, in accordance with the correspondence principle as introduced by Niels Bohr 
in 1923. We note furthermore that insisting that t' < t ensures that account is taken 
of retardation, corresponding to the Lifshitz definition of the retarded photonic Green's 
function (e.g. p2] §75) which is for t < t'. 

We make use of the Fluctuation-Dissipation theorem at zero temperature, rendered 
conveniently as 

i(E i {r)E k (r')) u = -3m{r !t (r,r»} (4a) 
i(^(r)i/ fe (r / )). = --^Curl^-Curl^m {T^r, r'; u)} , (46) 

Hq flfi ui 

|| Compared to Lifshitz et.al. T> = —hc 2 T/u! 2 , which is only a matter of definition. 



Casimir force on slab in cavity 



£ l^l £ g A ^2 £ g ,^ g £ 1,^1 




Figure 1. The five-zone geometry of a slab in a cavity. We have chosen z = at the 
left hand cavity wall. 



with the notation Curljfc = tijkdj [tijk being the Levi-Civita symbol and summation 
over identical indices is implied), Curl^ = e^d'^ where dj is differentiation with respect 
to component j of r', and fi' = fi{r'). The brackets denote the mean value with respect 
to fluctuations. The Casimir pressure acting on some surface is now given by the 
zz-component of the Abraham-Minkowski stress tensor, found by simple insertion to 
become [T71 [2 ljll 

(5) 

where a standard frequency rotation u = i( has been performed and the convenient 
quantities T E and T H have been defined according to 

T ?k( r i r '; w ) = F ik (r, r'; u); (6a) 
Ff k (r, r'- u) = ^Cm\aCm\' km T lm (r, r'; u). (6b) 

In (JSJ) only the homogeneous (geometry dependent) solution of (T2]) is included; the 
inhomogeneous solution pertaining to the delta function represents the solution inside 
a homogeneous medium filling all of space. This term is geometry independent, and 
cannot contribute to any physically observable quantity. Importantly, however, any 
such simplification from the full Green's function to its "effective" counterpart must 
only be made subsequent to all other calculations. 

The system is symmetrical with respect to translation and rotation in the xy-plane 
and we transform the Green's function once more: 

? (r,r';w) = J |^ e *^-*i) I (z,J;k ± ,u). (7) 

% The expression is generalised compared to the original reference to allow fi ^ 1. 
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Here and henceforth, the subscript _L refers to a direction in the rry-plane. In the k±,u 
Fourier domain one finds [SI [52] that the component equations (J2J) combine to (among 
others) the equations 

k 2 

{d 2 z - K 2 )g xx {z, z'\ fc_L, u) = —8(z - z) (8) 

and 

,2 



( d z - n 2 )g yy (z,z';k ± ,uj) 



5(z - z'), 



(9) 



which readily give us these two components in each homogeneous zone. We have defined 
the quantity k = {k\ — efito 2 /c 2 ) 1 / 2 . The final diagonal component is found by means 
of the relations 

ik± 



g zz (z, z';k±,u) 



g zx (z,z';kj 
g xz (z, z';kj 



K 

ik 



K z,z ; k^uj) + —^-^(^ 



:iOa) 



-^-d z g xx (z,z'; k ± ,u) 



9zx(z', z\ -fej 



(106) 
(10c) 



of which the first two are components of (T5]) and the last was shown by Lifshitz (e.g. 

inn- 

We return to the geometry of figure [TJ An important point to emphasize is that 
unlike certain authors in the past (e.g. [53]) we make no principal difference between 
the walls of the cavity and the slab; they are both made of real materials with finite 
permittivity and conductivity at all frequencies as is the case in any real experimental 
setting. The net force density per unit transverse area acting on the slab is found by 
first placing the source (i.e. z') in one of the gaps and calculate the resulting Green's 
function in this gap. This yields the attraction the stack of layers to the left and right 
of this gap exert upon each other. The procedure is then repeated with respect to the 
other gap region and the net force on the slab found as the difference between the two. 

The solution of ([S]) may be written down directly, yielding in the case where z' lies 
in the gap region to the left of the slab 

Ae KlZ z < 

Cie K ^ + C 2 e~^ z + Ge- K ^ z - Z '\ < z < a + 



E x e K2Z + E 2 e- K2Z 
Die KgZ + D 2 e~ KgZ 



< z < a + + b 
' + b < z < c 
z > c 



with G = —Kg/(2e g fj, g ). The "constants" A to E are ^'-dependent. From standard 
conditions of EM field continuity and fj¥q[6) one may show that g xx and (e/ K 2 )d z g xx are 
continuous across interfaces, giving a total of 8 equations which are solved with respect 
to C\ and C 2 yielding after lengthy but straightforward calculation (for details, cf. [54J) 
the solution in the left hand gap (0 < z < a + , denoted with superscript +) 

1 



9.i 



d 



TM 



2 cosh K g (z 



Kg{z + Z') 



A™ 



VM e -K g (z+z') 



+A1 1 



TM. 



-Kg(z + Z') 
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Here and henceforth the inhomogeneous \z — z'|-term has been omitted subsequent to 
other calculation as argued above. Foreknowingly, we have defined the key quantities 
1 



" q q ^ 



q = {TE, TM} 



(12) 



where 



± 



A 1? A 2? (1 
1 



A lq A 2q e~ 2 ^ a± \ 



A lq A 2q e- 2 ^ a± - A 2q (A 2 



q - — iq— zq- ^2q\^2q ~ ^lq 

using the single-interface Fresnel reflection coefficients 



A lo e- 2K9a± )e 



-2K 2 b 



A,, 



™i l/i,q 6 



l.q 



Ha 




1,2. 



(13) 

Note already how the quantity (d^)" 1 is a generalisation of the quantity d~ l as it was 
defined for the three-layer system by Schwinger et.al. [H [52] (dubbed A in the Lifshitz 
et.al. literature). In the limit k 2 — > k 9 we immediately get {d^)^ 1 — > (A lq 
i.e. the three-layer standard result for a cavity of width c with no slab. 
Following the above described procedure we get 
1 



2 e 2 Kg c _ -p 



2n g e g 



d 



TM 



2 cosh K g (z — z') — 



e Kg(z + z') 



A? 



TM 



-A! 



TM e -K s (z+z') 



+ A X 



TM e -K g {z+z') 



Exactly the same procedure as for g xx is followed to obtain the ^-component. One 
finds that g yy and fi~ l d z g yy are continuous across boundaries, giving 8 new equations 
solved as above to yield 



9yy(+) 



JXg_UJ_ 



1 



d 



TE 



D Kg(z + z') 



2 COsh Kg (z 



z - 



Af 



TE 



-A( 



TE -Kg (z+z') 



Af 



TE, 



-Kg(Z+Z f ) 



The results for the right hand (— ) gap is found by transforming the above results 
according to a ± — > a T and z — >• c — z. 

To obtain the force density on each side of the slab, the solutions are now inserted 
into One may show [51] that the terms depending on z + z' do not contribute to 
the force density (this is a subtle point which will be discussed further below). Upon 
omitting these terms, the right hand expressions are simply given by swapping + and 
— indices everywhere and we are left with the effective Green's function solution in the 
w, k i -domain: 



9xx(±) 
9yy(±) 
9xx(±) 



•— — j— cosh kJz — z') 

€n dn 



TM 

^-±-COShK a (z 



4 TE 
1 



Kg€g d^ M 



cosh K g (z 



'Ma) 



(146) 



;i4 C ) 



Upon insertion into (jSJ) we find the force on either side of the slab yielding the net 
Casimir pressure acting on the slab towards the right as 

TM 



J 7 °(a + ,a ;b,c) 



2^ 



dfc 



<j=TE 



1 1 

dz~d+ 



(15) 
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Naturally, the force will always point away from the centre position. Superscript 
here denotes that the expression is taken at zero temperature. The finite temperature 
expression, as is well known, is found by replacing the frequency integral by a sum over 
Matsubara frequencies according to the transition 



ft / ^f(K) - k B TJ2 f(Km); Km = i(2nk B T/H) 

J{i m=0 



rn 



yielding 



U rji 00 ■ f00 ™ / 1 1 \ 

11 m=0 J ° q=TE \ % % ' 

The prime on the summation mark denotes that the zeroth term is given half weight as 
is conventional. 

Rather than painstakingly solving the eight continuity equations to obtain the 
Green's function as above, the result (fT5l) is found much more readily using a powerful 
procedure following Tomas as presented recently by one of us [20J . The above result was 
obtained by Tomas [55] presumably using this procedure. It was worth going through 
the above calculations, however, for the sake of shedding light on some in our opinion 
non-trivial details which are often tacitly bypassed. 



3. Casimir measurement by means of an oscillating slab 

Equation (fT6l) may be written on a more handy form in terms of the distance 5 from 
the centre of the slab to the midline of the cavity. We introduce the system parameter 
h = c—b = a + +a~ and substitute according to a = h/2±5. With some straightforward 
manipulation we are able to write ffT6l) as 

k R T^' r „ , ^ A,sinh2/c„5 



with 



m=0 Jl) q=TE 11 9 



A g = 2A lq A 2q (l-e- 2K2b )e- K *\ 

B q = l- Al q e~ 2K2b + A 2 lq (A 2 2q - e -^ e -2 Ka h 

We write the force on the slab at finite temperatures as a Taylor expansion to first 
order in 5 as 

F T (8; b, c) = ai 6 + 0(8 3 ) (18) 

with 

a 1 (T;b,c) = 2 -^Yl / d ^-^E^V" ( 19 ) 

71 m=0 J q =TE^1 A 1 
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Figure 2. The slab oscillates about the cavity midline. We imagine a spring is 
attached to the slab exercising a Hooke-force towards the equilibrium position. 



Assume now the slab is attached to a spring with spring constant k per unit 
transverse area. For small 5 we may assume the slab to oscillate in a harmonic fashion 
(assuming k > a\ now) with frequency given by Newton's second law as 



Q = Q - AQ(T) 



k — ai(T) 



'in 



where f2 = \Jk/m and m is the mass of the slab per unit transverse area. In the case 
that > ai we get 

Afi(T)^^2 = fi n ai(T) 



2\/km " 2k 

We show by numerical calculation in chapter [5] how the Taylor coefficient a\ (T) varies 
significantly with T rendering an oscillating slab-in-cavity setup possibly suitable for 
future experimental investigation of the true temperature dependence of the Casimir 
force. 

The setup as described is somewhat reminiscent of the setup currently employed 
by Onofrio and co-workers in Grenoble [50] where plates mounted on a double torsion 
balance are attracted to a pair of fixed plates. In their planned experiment, the distance 
from plate to wall will however kept constant during force measurements. Indeed, a 
double torsion balance might be one way of envisioning an experimental realisation 
essentially equivalent to the system described (if thickness corrections are neglected) 
if the plates are mounted such that when one pair of plates approach each other, 
separation is increased between the pair on the opposite side of the pendulum. An 
even closer relative might be the recent experiments in Colorado where perturbations 
of the eigenfrequency of a magnetically trapped Bose-Einstein substrate in the vicinity 
of a surface provides a sensitive force measurement technique [5] [301 EE]- Both of these 
experiments involve a plate (in the widest sense) attracted to a wall on only one side; 
an "open cavity". 

While a one-sided configuration is possibly experimentally simpler, there are two 
physical advantages of the sandwich geometry as presented here: the frequency shift 
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Afi(T) is essentially twice as large using a closed cavity and, perhaps more importantly, 
in a symmetrical geometry the harmonical approximation (JF T oc 5) is accurate for larger 
deviations 8 from the equilibrium position than is the case for an open geometry. These 



points are elaborated further in Appendix A 



4. Fundamental discussion: two-point functions and Green's functions 

In the standard Casimir literature there are two famous and somewhat different 
derivations of the classical Lifshitz expressionEl, namely that of Lifshitz and co-workers 
in 1956-61 [T5l [TT] and that of Schwinger and co-workers some years later [8j [52]. The 
two both make use of a Green's two-point function but in two different ways which upon 
comparison seem somewhat contradictory at first glance. Understanding how they relate 
to each other is not trivial in our opinion. 

In order to calculate the force acting on an interface between two different media, 
both schools calculate what in our coordinates is the zz component of the Abraham- 
Minkowski energy momentum tensor as described above using the Green's function 
through the fluctuation-dissipation theorem as in (j4"a] b). Lifshitz argues as recited above 
that in his formalism some terms of the Green's function (those dependent on z + z') 
make no contribution to the forced. These are consequently omitted, leaving an effective 
Green's function. Schwinger et.al., however, make use of the entire Green's function 
ultimately arriving at an expression similar to (jSJ) in which the z + z' terms are included 
and indeed necessary in order to reproduce Lifshitz' result. The \z— z'\ dependent source 
term is geometry independent and eventually omitted in both references. 

To solve the paradox we recognise one important difference between the two 
procedures: Lifshitz takes the limit r — > r' so that r and r' are both on the same 
side of one of the sharp interfaces, whereas in Schwinger's method, r is on one side 
whilst r' is on the other. By using continuity conditions for the EM field, calculations 
can be carried out with analytic knowledge of the Green's function only on one side of 
the interface in both cases, thus masking this principal difference. Remembering that 
T zz is the density of momentum flux in the z-direction, the physical difference between 
the methods is that whilst Lifshitz calculates the force density as the net stream of 
momentum into one side of the interface, Schwinger et.al. 's expression represents the 
entire stream into one side minus the entire stream out of the other side. Due to 
conservation of momentum, the procedures are physically equivalent. 

The question remains how to interpret the terms dependent on z+z'. Arguably, the 
absolute value of such terms must be arbitrary, since they will depend on the position of 
an arbitrarily placed original Furthermore, since these terms cancel each other perfectly 

+ By "Lifshitz force" is henceforth meant the Casimir force between two plane parallel 
(magneto)dielectric half-spaces separated by a medium different from both. By the "Lifshitz expression" 
is meant the mathematical expression for this force as derived by Lifshitz and co-workers [21 . 
* This is shown formally in [54] . 

jj The notion of arbitrarily large energy densities, of course, is not foreign to Casimir calculations; 
Casimir's original calculation involved the difference between the apparently infinite energy density of 
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Figure 3. Contributions to g in a gap between two bodies in the optical visualisation. 
The distance between the bodies is a. Each term has a weight factor as shown on the 
right hand side. The sum of the infinitely many reflections of a g-polarised wave is 

in (jSJ), one may think of them as representing an isotropic flux of photonic momentum, 
flowing in equal amounts in both directions along the z-axis, giving rise to no measurable 
effect inside a homogeneous medium. 

Schwinger, however, insists r and r' lie infinitesimally close to either side of an 
interface. While the z + z' terms cancel each other when all calculated in the same 
medium, their values depend on e and fi, so when e ^ e' or fi ^ fi', their net contribution 
is finite. 

This is exactly made up for in Lifshitz' approach by the fact that a sudden change in 
permittivity and permeability (such as at an interface between a dilute and an opaque 
medium) causes some of the radiation to be reflected off the interface in accordance 
with Fresnel's theory. Thus although z and z' both lie inside the same medium, there 
is a net flow of momentum either out of (attractive) or into (repulsive) the gap giving 
rise to a Casimir force. Such an analysis of the use of Green's functions gives way 
for an understanding of how three different representation of the Casimir effect come 
together; the derivation by Lifshitz starting from photonic propagators in quantum 
electrodynamics, that by Schwinger et.al. based on Green's function calculations from 
classical electrodynamics and a third approach based on Fresnel theory which we may 
refer to as the "optical approach" (originally in form of non-retarded Van der Waals 
theory (56J [57], recently revisited by Scardicchio and Jaffe, see [58J and references 
therein) . 

We showed that the factors (c^) _1 were generalised versions of the factors denoted 
the zero-point photon field in the absence and presence of perfectly conducting interfaces. 
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d 1 and {d!) 1 in Schwinger et.al.'s theory for the three-layer model. These are both 
special cases of a more general quantity 

d q 1 - r qh r qR e- 2K s a 

pertaining to a gap of width a separating planar bodies to the left (L) and right (R) 
of it whose Fresnel reflection coefficients are r g L and r qR respectively. If the media 
are infinitely large and homogeneous media indexed 1 and 2 respectively, say, r q ^ and 
r qR are simply — A lq and — A 2q from f|T3|) ; if the bodies are more complex, e.g. has a 
multilayered structure, their corresponding Fresnel coefficients will be more complicated. 
This is discussed in detail in [20]. An EM plane wave with momentum Kk is described 
as e 1 ( kj -' r±+fcz2; ). In medium g, furthermore, k z = in g according to Maxwell's equations, 
i.e. the wave is evanescent in the ^-direction if k\ > e g fi g uj 2 /c 2 (otherwise propagating). 
After frequency rotation u 2 — > — C, 2 this is always true (k± is assumed real), so every 
wave is described as an evanescent wave. The attenuation of an EM field of frequency 
i( propagating a distance I along the z-axis in medium g is exp(— n g l), so one readily 
shows that d~ x is the sum of relative amplitudes of the electric fields having travelled 
all paths starting and ending at the same ^-coordinate and with the same direction: 

1 oo 

T = r qh r qR e- 2 ^ a + (r^Re" 2 ^) 2 + ... = ^V^^T- 

H n=l 

An illustration of this is found in figure [31 Since the phase shift from propagation in 
the _L direction is disregarded in this respect, one might think of d~ l as a sum over all 
closed paths, parallel to the 2-axis and starting and ending in the same point. 

Considering again the expressions for the complete Green's functions g X x,9yy an d 
g zz in section 2, we see that the last term of all three components are the only ones 
not multiplied by a factor d~ x (indices ± suppressed). Since this factor is the only 
part of g containing geometry information, the last term is geometry independent, and 
can obviously make no contribution to a physical force. Hence: all contributing terms 
are proportional with d~ x which leads us to the conclusion that the Casimir attraction 
between bodies on either side of a gap region at a given temperature depends solely on 
the extent to which some EM field originating in the gap, stays in the gap. 

To sum it all up, we argued that Schwinger's classical Green's function as introduced 
is the exact macroscopic analogy of Lifshitz' QED propagator according to Bohr's 
correspondence principle. In its Fourier transformed form it expresses the probability 
amplitude that an electric field which has transverse momentum ^.k^, energy ftio and 
coordinate z' will give rise to a field of the same energy and momentum at z. When then 
z and z' are only infinitesimally different, the only ways this can happen by classical 
reasoning is that the two are in fact exactly the same (corresponding to the \z — z'\- 
dependent source term) or that the field has been reflected off both walls once or more. 
This is what figure [3] demonstrates. 
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Force on centre slab in 5-zone system 

10 2 F 1 1 1 1 




10 1 1 1 1 1 1 1 

200 400 600 800 1000 1200 

Distance 8 from centre line/nm 



Figure 4. The force on an Al slab in a vacuum-filled cavity between Al walls. 5 is the 
distance from the centre of the slab to the midline of the cavity. For negative 8 one 
gets the antisymmetrical extension of the graph. 

5. Numerical investigation and temperature effects 

For our numerical calculations, we have used permittivity data for aluminium, gold 
and copper supplied by Astrid Lambrecht (personal communication). For ease of 
comparison, aluminium is used in figures throughout; all variations aqcuired by replacing 
one metal by another are of a quantitative, not qualitative nature, and are not included 
here. In all our numerical investigations, we have assumed non-magnetic media, i.e. 

Ml = ^2 = l^g = 1- 

As an example of a dielectric, we have chosen teflon fluorinated ethylene propylene 
(teflon FEP) because its chemical and physical properties are remarkably invariant with 
respect to temperature. Permittivity data for teflon FEP are taken from [59] . 

Figure H] shows the Casimir force acting on a relatively thick aluminium slab in a 
cavity as a function of 5. For negative values of S the situation is identical but the force 
has the opposite direction. We have chosen a gap width of 3/xm and a slab thickness 
of 500nm. These values are not arbitrary: First, the relative temperature corrections 
of the Casimir force are predicted to be large at plate separations of 1-3/mi, so a slab- 
to- wall distance in this region is desirable (here h/2 = 1250nm). Secondly, choosing 
the slab significantly thicker than the penetration depth of the EM field makes the five- 
zone geometry instantly comparable to the well-known three-zone Lifshitz geometry of 
two half-spaces; for slabs of a good metal thicker than ~ 50nm there is virtually no 
difference between the five-zone expression as derived above and that which one would 
acquire applying the standard Lifshitz expression to each gap in turn and finding the 
net force density on the slab as the difference between the two. 
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Force on centre slab in 5-zone system 
relative to ideal result 
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Figure 5. The force on an Al slab in a vacuum-filled cavity between Al walls relative 
to its value for ideally conducting slab and walls, equation (|20|) . S is the distance 
from the centre of the slab to the midline of the cavity. For negative 6 one gets the 
symmetrical extension of the graph. 



Figure [5] shows the net vacuum pressure acting on the slab relative to Casimir's 
result for ideal conductors, 



(h/2 - S) 4 [h/2 + 5) A \ ' v ' 

In such a plot we see clearly how a slab and cavity set-up might be suitable for 
measurements of temperature effects; whereas such effects are small for very small 
separations, they grow most considerable near the centre position where slab-to-wall 
distance is in the order of a micrometre. 

An altogether different result is obtained upon replacing metal with a dielectric in 
both walls and slab. In figure M the same calculation as in figure has been performed 
with both slab and walls of teflon FEP. Casimir experiments using dielectrics were 
proposed by Torgerson and Lamoreaux [60] where the use of diamond was suggested. 

It is important to note here that we have not taken into account variations of the 
dielectric properties of teflon FEP with temperature; much as teflon FEP is renowned 
for its constancy in electrical and chemical properties over a large temperature range and 
is used in space technology for this very reason, one must assume there are corrections at 
extremely low temperatures. We shall not enter into a discussion on material properties 
here; the point to take on board is rather that temperature effects are found to be very 
large indeed near the centre position, a fact that does not change should the calculated 
values be several percent off. This strongly indicates that the use of dielectrics in Casimir 
experiments could be an excellent means of measuring the still controversial temperature 
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Force on centre slab in dielectric system 
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Figure 6. The force on a teflon FEP slab in a vacuum- filled cavity between teflon 
FEP walls relative to the result for ideally conducting slab and walls, equation 
Note that dielectric properties are assumed constant with temperature. 



dependence of the force. 

We note furthermore that whilst for metals the force decreases with rising 
temperatures, the opposite is the case for the dielectric. Mathematically this is readily 
explained from e.g. (1161) . Temperature enters into the expression in two ways; first, 
each term of the Matsubara sum has a prefactor T, secondly the distancing of the 
discrete imaginary frequencies increases linearly with T. The first dependence tends to 
increase the force with respect to T whilst the other decreases it (bearing in mind that 
the integrand, which is proportional with exp(— K g h), decreases rapidly with respect to 
£ for ( larger than roughly the m = 1 Matsubara frequency). As temperature rises, 
thus, the higher order terms of the sum quickly become negligible, leaving the first few 
terms to dominat© In the high temperature limit, m = becomes the sole significant 
term and the force becomes proportional^ to T. This is true for metals and dielectrics 
alike, but while the trend is seen at low temperatures for dielectrics, for metals the 
T-linear trend typically becomes visible only at temperatures much higher than room 
temperature. In metals the low (nonzero) frequency terms are boosted since e« ^> e g 
for ( much smaller than the plasma frequency, in which case reflection coefficients |Aj g | 
approximately equal unity. The first few Matsubara terms thus remain significant as 
temperature rises, countering the T-proportionality effect, at the same time as each 
m > term decreases in value as the Matsubara frequencies take higher values, allowing 

ff The same phenomenon for increasing distances rather than temperatures is treated in |62j . 

f For the three-layer Lifshitz set-up, the zero term and thus the force becomes proportional to T/a 3 

where a is the gap width, as shown formally in [llj . 
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Force density on slab in a cavity 



Cavity width: 3nm 
Slab thickness: 500nm 
Slab/cavity material: 
Gap material: vai 
Temperature 300K 
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Figure 7. The Casimir force density on the slab of figure Q] as a function of the 
distance 5 from the centre of the slab to the cavity midline as compared to its first 
order Taylor expansion, equation (|18p at temperature 300K. 



the resulting force to decrease with increasing temperature. 

Figure [7] shows the force acting on the slab in the previously described geometry 
(such as plotted in figure HI) as well the first order Taylor expansion. The figure gives 
a rough idea as to the size of the central cavity region in which one may regard the 
force density as linear with respect to 5. With the system parameters as chosen we see 
that, depending on precision one may allow oscillation amplitudes 6 of several tens of 
nanometres, a length which is not small relative to the system. 

The first order Taylor coefficient itself has been calculated and plotted in figure M 
for aluminium and teflon FEP slabs in an aluminium cavity. These are furthermore 
compared to Casimir's ideal result ( 1201) whose first order Taylor coefficient is readily 
found to be 

lQhcir 2 . 
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3.3283 • 10" 25 Nm 2 • h' 



(21) 



6. The effect of finite slab thickness 



As measurements of the Casimir force have become drastically more accurate over the 
last few years, with researchers claiming to reproduce theoretical results to within 
1% 0, [22], it is well worth asking whether any theoretical calculation may rightly 
claim such an accuracy. A point of particular interest in this respect is the strong 
dependence of the Casimir force on the permittivity of the media involved. The 
permittivity data for aluminium, copper and gold supplied by Lambrecht and Reynaud 
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Taylor expansion of cavity-slab force 



Aluminium slab 
Teflon FEP slab 
Ideal conductors 



Cavity material: Aluminium 
Gap material: Vacuum 
Cavity width: 3um 
Slab thickness: 500nm 
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Figure 8. The first order Taylor coefficient of equation (| 18[) for aluminium and teflon 
FEP slabs in an Al cavity. The horizontal dotted line is the coefficient pertaining to the 
Casimir result for ideal conductors (both slab and walls), equation (f2Tj) . One should 
note that the dielectric properties of the materials at extremely low temperatures are 
not known. 



were calculated by using experimental values for the susceptibility at a wide range of 
real frequencies (approx. 1.5 • 10 14 rad/s< lu < 1.5 ■ 10 19 rad/s), extrapolating towards 
zero frequency by means of the Drude relation (for small u < approx 1.5 • 10 14 rad/s). 
e{uj) was subsequently mapped onto the imaginary frequency axis invoking Kramers- 
Kronig relations numerically. Thus, although matching theoretical values (Drude mode 
1) excellently for imaginary frequencies up to about 10 15 rad/s [21], the data have 
intrinsic uncertainties. Recently, Lambrecht and co-workers addressed the question 
of the uncertainty related to calculation of the Casimir force due to uncertainty in 
the Drude parameters used for extrapolation, found to add up to as much as 5%, 
considerably more than the accuracy claimed for the best experiments to date [61J. 

The effect of the "leakage" of vacuum radiation from one gap region to the other in 
our 5-zone geometry is worth a brief investigation in this context. Excepting the zero 
frequency term, it is unambiguous from e.g. the definition of k that when the slab is 
metallic, the factor exp(— 2k 2 6) is small compared to unity for sufficiently large values 
of b, due to the large values of e 2 for all important frequencies ic§ Let us regard one 
of the gap regions, of index ±. With some manipulation one may expand ( 1121) to first 

| The m = term is a subtle matter we shall not enter into here. For a recent review, see e.g. [29] and 
references therein. 
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order in the factor exp(— 2/t 2 6) to find the pertaining quantity 
J_ A l9 A 2g e- 2 ^ a± 

-2 K2b A lg (l-Al g )e-^ a± A 2q - A lq e~^ 

(1 - A lq A 2q e-^y ' 1 - A lq A 2q e~^ + ^ j ' ^ 
The first term is immediately recognised as giving the Lifshitz expression for the Casimir 
attraction between two half-spaces of materials 1 and 2 separated by a gap of width 
and material g, and the second term is the first order correction due to penetration of 
radiation through the slab. 

In terms of 5 we may write in the case where exp(— 2n 2 b) <C 1 for all relevant 
frequencies (again subsequent to some manipulation) the force on the slab as .F T (<5) ~ 
J 7 ^ + AjF T where 

k B T^' [°° „ , ^ A qh smh2K g S 



B t — A oL cosh 2n a S 

m=0 " u <j=TE 9ij 9 

is the result using the Lifshitz expression on both gaps and taking the difference; here 
A q L = 2A lq A 2q e~ Kgh and B qL = 1 + A\ q A\ q e~ 2Kgh ; (23) 



and 



. rp oo ' oo TM 

A^ T (<5;M) = -^ij] / dA; ± A; ± «: g ^e- 2K2fe x 

™ — n ^ j-tt 



m=0 ^ u <7=TE 

x A qL (B qL - A\ q - A\ q e-^ h ) sinh2^ 

- -4gL COsh 2k 9 £) 2 

The factor exp(— 2k 2 6) and consequently the first order correction, is very sensitive 
with respect to even small changes in e 2 (i£). For very thin slabs (b < 50nm) and small 
cavities, the correction could be in the order of magnitude of the currently claimed 
measurement accuracy. Furthermore, we see that the integrand of (|2"3|) depends on e 2 
in an exponential way. In conclusion: to the extent that the thickness correction is of 
significance in an experimental measurement, exact knowledge of the permittivity as 
a function of imaginary frequency is of the essence. In such a scenario, approximate 
knowledge of the dispersion function could effectively limit our ability to even calculate 
the force with the precision that recent experiments claim to reproduce theory [22J. 

In the case of dielectrics, the correction is almost two orders of magnitude larger 
and should be readily measurable. Experiments in a geometry involving dielectric plates 
of finite thickness might even be a possible means of evaluating the correctness of the 
dielectric function employed. 



7. Conclusion and final remarks 



The main conclusion from the work presented is that from a theoretical point of view 
the five-zone setup (figured!) as discussed could be ideal for detection of the temperature 
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Finite slab thickness correction 




Figure 9. Thickness correction due to "leakage" of radiation through a thin slab in 
a small cavity. The absolute value of the correction is approximately exponentially 
decreasing with the thickness b. 



Finite slab thickness correction, teflon FEP 
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Figure 10. Thickness correction due for teflon FEP. 
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dependence of the Casimir force when the wall-to-slab distance is in the order of 1/im. 
One method as suggested is a measurement of the difference in the eigenfrequency of an 
oscillating slab in the absence and presence of a cavity. 

When metal is replaced by a dielectric in slab and walls, relative temperature 
corrections become much larger, suggesting that using dielectrics whose dielectric 
properties vary little with respect to temperature be excellent for such measurements. 

Our treatment of the effect of finite slab thickness shows that the effect of finite 
thickness varies dramatically with respect to the properties of the materials involved, 
specifically e and //. Much as the effect is generally quite small for metals, to the extent 
such effects do play a role even a moderately good estimate of their exact magnitude 
requires very accurate dielectricity data for the material in question. This is but one 
example of the more general point that the still considerable uncertainties associated 
with the best available permittivity data for real materials call for soberness in any 
assessment of our ability to numerically calculate Casimir forces with great precision. 

Finally, a couple of remarks: it ought to be pointed out that our proposed method of 
investigating the thermal Casimir force via observing the oscillations of a slab in a cavity, 
can be classified as belonging to the subfield usually called the "dynamic Casimir effect" . 
The use of mechanical microlevers has turned out to be very effective components for 
high sensitivity position measurements, of interest even in the context of gravitational 
waves detection. As for the basic principles of the method see, for instance Jaekel et 
al. [63] with further references therein, especially Ref. [M]. For more recent papers on 
microlevers, see Refs. [65l 1661 l6? t 168]. 

Moreover, we note the connection between our approach and the statistical 
mechanical approach of Buenzli and Martin [69|. These authors computed the force 
between two quantum plasma slabs within the framework of non-relativistic quantum 
electrodynamics including quantum and thermal fluctuations of both matter and field. 
It was found that the difference in the predictions for the temperature dependence of 
the Casimir effect are satisfactorily explained by taking into account the fluctuations 
inside the material. Their predictions for the force are in agreement with ours. 

Appendix A. Closed geometry versus open configuration 

We will demonstrate briefly the two physical properties favouring a closed cavity 
configuration (figure [1]) as compared to an open configuration in which a similarly 
oscillating plate is held in equilibrium by an external spring system. For the purpose of 
comparison we will disregard effects due to finite plate thickness, so that the net force 
experienced by a slab in a cavity is the difference between standard Lifshitz forces on 
both sides, whist that between plate and wall in a one-sided geometry (like figure [H but 
with the right hand wall removed) is simply the Lifshitz force. For separations of some 
hundred nanometres or more, the Lifshitz force varies as J-"h{d) oc d~ 4 . Thus the net 
force on the slab in the sandwich geometry attached to a spring of spring contstant k 
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Figure Al. The relative correction to the linear Taylor expansion, J-'uu, of the Casimir 
force near equilibrium position for the open and closed geometries plotted for positive S. 
Calculations assume aluminium plate and walls, equilibrium plate-to-wall separation of 
a =1250nm in both configurations and temperature 300K. For the sandwich geometry, 
T is given by ([16]) whilst in the open geometry T ~ T\,{a + 8) — T\,{a) with the 
standard Lifshitz expression for the attraction between two half-spaces. 



per unit transverse area is (we assume k > ^{a) /5 as before) 
Sandwich = Fi{a + 5)- F h (a - 8) - k8 

_ n*) m h 

{1 + 5/ay (1-5/a) 4 

= -[ka-8 |^ L (a)|] - + 40 |^ L (a)| + ... (A.l) 
a a 6 

where a = h/2 is here the distance from slab to wall in equilibrium position (J r h{d) < 0). 

Now consider an open configuration in which a plate is held in equilibrium by an 
external spring, also of spring constant k per unit transverse area. Assume that the 
forces are in equilibrium when the plate is a distance a from the cavity wall. The net 
force on the plate is 

^opon = ^L(a + S) - ^l(o) - k5 

5 <5 2 

= - [ha- A\T h {a)\}- - 10|JF L (a)|— + ... (A.2) 

There are thus two properties that favour the closed geometry. First, the first- 
order perturbation of the spring constant is twice as large and second, that the leading 
order correction to the harmonical approximation (JF T oc S/a) is cubical whist it is 
quadratic for the open configuration^. The closed geometry thus allows considerably 

§ Note that the specific assumption T^d) cx dr A is not necessary for either of these results; they 
pertain almost exclusively to geometry. A more general power J-h(d) oc d~ a , a > 0, say, gives the same 
properties. 
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larger deviations from equilibrium position at a given accuracy without taking non- 
harmonic effects into account. In figure |A"T1 this is demonstrated by plotting the relative 
nonharmonic correction as a function of 5 at a separation 1250nm. In accordance with 
our results, the relative correction {T — T\j^)jJ- is approximately linear for an open 
geometry (« — §-) and approximately quadratic for a sandwich (« 5^). 
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